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SYSTEM AND METHO " FOR CHARACTER^ A SAMPLE BY LOW- 
FRFm IFNCY SPECTRA 



CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is a continuation-in-part, and claims the benefit of, U.S. 
Patent Application No. 10/112,927, filed March 29, 2002, entitled APPARATUS 
AND METHOD FOR MEASURING MOLECULAR ELECTROMAGNETIC 
SIGNALS WITH A SQUID DEVICE AND STOCHASTIC RESONANCE TO 
MEASURE LOW-THRESHOLD SIGNALS (Attorney Docket No. 
38547.8007.US00), International Patent Application No. PCT/US03/09544, filed 

u „ a onno «„««arf RYSTFM AND METHOD FOR CHARACTERIZING A 

SAMPLE BY LOW-FREQUENCY SPECTRA (Attorney Docket No. 

Mill IMAAIiiA^i 

38547 8005 WO00), and International Patent Application No. rui/uouon ioa-, 
filed April 18, 2003, entitled SYSTEM AND METHOD FOR SAMPLE DETECTION 
BASED ON LOW-FREQUENCY SPECTRAL COMPONENTS (Attorney Docket 
No. 38547.8006.WO00), all incorporated herein by reference. 

BACKGROUND 

There are a variety of spectroscopic tools for characterizing atomic or 
molecular compound. These include, but are not limited to, x-ray, UV, visible-light, 
infrared and microwave spectroscopy, and nuclear and electron spin resonance 
(NMR and ESR) spectroscopy. In general, spectroscopic tools are useful for at 
least four drfferent type of chemical-analytical problems: first, to characterize an 
atomic and molecular compound according to its spectrographs features, e.g., 
spectral components; second, to determine the atomic composition of a 
compound, according to the spectral characteristics of atoms making up the 
compound; third, to determine 2-D or 3-D conformation of a molecular compound 
according to the spectral characteristic of atom-atom interactions in the compound; 
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**** . ^ tools provide some unique advantages) in 

Most existing spectroscopy tools pm ^ ^ 

ttms ot sens* ^^2^Zl -* * is gene* 

pertinentspectroscopictoolsaspossible. 
SUMMARY 

a «n P ct a method of characterizing spectral 
, The invent includes ,n onea a ^ — ^ fc 

e m *s to n features ot a sample ^^ liBMlll rt*-<«. 
secular motion within the sample. The metho ^ ^ ^ ^ ^ 

.-~, 3l uhere F* T is the total sample count S, PP 
*"'•*" . ,^ p » a i Fast Founer Transform of the time oom 

^domain resoluhon tola Real Fast Fo ^ ^ ^ ^ ^ ^ ^ 

s ,„a, sanded at sampli* r**J ^ J ^ 8]gnal ^ performs a 
p^ms^S/nsamp^rom ^ ^ ^ . ^ 

Real Fast Fourier Transform (RFFT) on tn v ^ forthe 
sigoal is calculated from the normalized ^ ' ^ ^ ^ e measure d power a. 

a „d c is chosen such that the *>ta ^ ^ ^ 

« 20-50% of the — * „ each even, Pin 

ste psarerepea W n«n.s.genera.in,ah^m«* n 

fram the RFFT in f conespondrng-freguenoy p o ^ ^ f 

power bins by n, to y w,u 10/9/0 3 
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[0007] 



histogram, the average power in each bin. The method may further include 
idenWyfng those bins in the histogram that have an event count above a grven 

threshold and an average power. 

Also disclosed is a low-frequency spectral signature associated w»th a 
material of interest comprising a list of frequency components in the DC-50 kHz 
frequency range that are generated by the above method. The frequencies in the 
,ist may be identified from a histogram of the number of sample^lependent 
stochastic events occurring at each of a plurality of spectral increments withm a 
selected frequency range between DC and 50 kHz. 

In another aspect, the invention includes an apparatus for interrogating a 
sample that exhibits low-frequency molecular motion. The apparatus includes a 
magnetically and electromagnetically shielded container adapted for receivmg the 
sample, an adjustable-power source of Gaussian noise for injection into the 

Joiner and a detector for detecting an 

sample, wnn me sainyw «« 

aiectromaqnetic SmeOomain signal composed of sampte source radiatton 
superimposed on the injected Gaussian noise. An electronic computer rn ma 
apparatus receives the time-domain signal from the detector, and processes the 
signal to generate a spectral plot that displays, at a selected power setting of the 
Gaussian noise source, low+equency spectral components characteristic of the 
sample in a selected frequency range between DC and 50 kHz. 

In one general embodiment, the electronic computer includes a signal 
analyzer that functions to (i) calculate a series of Fourier specfra of the time- 
domain signal over each of a plurality of defined time periods, in a selected 
frequency range between DC and 50 kHz, and (ii) average the Fourier spectra. 
Typically, at least five Fourier spectra are calculated, each taken over a 1-5 

second time-domain interval. 
mm in another general embodiment, the electronic computer includes machine- 
readable code operable to cany out the method described above for generating a 

histogram of spectral events. 
m „ The source of Gaussian noise in me apparatus may be an adiustable-powe. 
Gaussian noise generator and a Helmholz coil whkth is contained within the 
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[0013] 



[0014] 



magneto cage and the Faraday cage, and which receives a selected noise output 
signal from the noise generator in me range 100 mV to 1 V. The injector ,s 
designed to inject Gaussian noise into the sample at a frequency, for example, 

between DC and 8 kHz. 

The detector in the apparatus may be a first-derivative superconducting 
gradiometer which outputs a current signal, and a SQUID operative* connected to 
the gradiometer to convert the current signal to an amplified voltage signal. 

The container in the apparatus may include an attenuation tube having a 
sample-holding region, a magnetic shielding cage surrounding the region, and a 
Faraday cage contained within the magnetic shielding cage and also surrounding 
the region. In this embodiment, the source of Gaussian noise may include a 
Gaussian noise generator and a Helmholz coil which is contained within the 
magnetic cage and the Faraday cage, and which receives a noise output signal 

.... „j fnrfhor induces, for use in removing 
from the noise genwaiui, a..u ........ 

stationarv noise components in the time-dependent signal, a signal inverter 
operative* connected to the said noise source and to said SQUID, for receiving 
Gaussian noise from the noise source and outputting into said SQUID, Gaussian 
noise in inverted form with respect to the Gaussian noise injected into the sample. 

In still another aspect, the invention includes a method for interrogating a 
sample that exhibits low-frequency molecular motion. In practicing the method, 
the sample is placed in a container having both magnetic and electromagnetic 
shielding, and Gaussian noise is injected into the sample at a selected noise 
amplrtude. An electromagnetic time-domain signal composed of sample source 
radiation superimposed on the injected Gaussian noise, is recorded, and from this, 
a spectral plot that contains, at a selected power setting of the Gaussian noise 
source, low-frequency, sample-dependent spectral components characteristic of 
the sample in a selected frequency range between DC and 50 kHz is generated. 
The steps are repeated at dffferent selected noise amplitudes until a plot showing 
a maximum or near maximum number of spectral components characteristic of the 
sample are generated. 
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[0019] 



,n one embodiment, the spectral plot is generated by (i) calculating a senes 
of Fourier spectra of the time-domain signal over each ofapluralrty of defined time 
periods, in a selected frequency range between DC and 50 kHz, and (ii) averaging 

the Fourier spectra. 

In another general embodiment, the spectral plot is generated by the 

histogram method above. 

These and other objects and features of the invention will become more 
yiy apparent when the following detailed description of the invention is read n 
conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is an isometric view of one embodiment of a molecular 
electromagnetic signaling detection apparatus formed in accordance wrth one 
embodiment of the present invention, 

n... 0 is an enlaraed, detail view of the faraday cage and its contents 

• ■» 

shown in Figure 1; and 
„ Ffcure 3 is an enlarged, cross sectoral view of one of .he attenuaflon tubes 

shown in Figures 1 and 2. 
„„ Figure 4 is a ooss-seotion view of the faraday cage and its contents shown 

in Figure 2. 

[0022] Figure 5 is a cross-section view of an alternative embod.ment of the 

invention shown in Figures 1 through 4. 
l0 o 231 Figure 6 is an enlarged, detail view of the frames supporting the co.ls of the 

Helmholtz transformer described herein. 
[0024] Figure 7 is a diagram of an alternative electromagnetic emission detection 

system. t 
{Ssm Figure 8 diagram of the processing unit included in the detection system of 

the above Figures. 

, 00261 Figure 9 is a diagram of an alternative processing unit to that of Frgure 8. 
m Figure 10 is a flow diagram of the sfcnal detection and pn.cess.ng 
performed by the present system. 
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m] Figures 12Aand 12B are spe |" ^| Fourier transforming a non- 
530 Hz, for a sample of saturated ^ '^^^ Courier transforming a cross- 
correlated contain sample sgnat (12A). and 
stated sampje spe*™ >m ^ ^ ^ 

^ T ir^ I su*. gene^d b, Fouder — a 
530 Hz, for a sample of alKyt ein transforming a 

non-correlated time^main sample agnal (13A), and Fou 

^ Ftgures 14A-14F are M P ^ ^ (14B) , . 

Hz, to samples of deumzed water 14AM ^ 
button of 1% NaCI in deionized water (14C), a satur 

... ..—-M4E). and nn sample (14FV 
aikyi ether »»»» »' — ' . enw ,„, (eqion between 500 and 

« 3 ° f 7,7^,4 , nrillion CM. 1:100 mlfon (15D), 
increasing w/v dMons «• ««» < 15B) ' B ^ generated 

v ,0 billion (15E and 15F), where me spe«ra , F^e « 



[0034] 
100351 



^LiTBisane^emenfofaporhonofFigurelTA 
1161 ^L.^seohonal^viewd.F.gureim 
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(0041] Figures 22A-22D are histogram spectra of a sample taken at four different 
noise power levels. 

[0042] Figures 23A-23C are computer screen shots displaying a user interface for 
generating and displaying a spectral plot histogram, with Figure 23C being in 
color. 

[0043] In the drawings, identical reference numbers identify identical or 
substantially similar elements or acts. To easily identify the disc 

J0044] Discussion of any particular element or art, the most significant digit or 
digits in a reference number refer to the Figure number in which that element is 
first introduced. 

DETAILED DESCRIPTION 



[0045] 
rnrufii 



I. Definitions 

. ii. . X-11...J — JoCnMnne ■ iniaee inrtiratari otherwise. 

I ne terms Deiow nave u re iunuwh iy v.^... .... 

"Sample that exhibits molecular rotation" refers to a sample material, which 
may be in gaseous, liquid or solid form (other than a solid metal) in which one or 
more of the molecular compounds or atomic ions making up or present in the 
sample exhibit rotation. 

"Magnetic shielding" refers to shielding that inhibits or prevents passage of 
magnetic flux as a result of the magnetic permeability of the shielding material. 
"Electromagnetic shielding" refers to, e.g., standard Faraday 

electromagnetic shielding. 

"Time-domain signal" or 'time-series signal" refers to a signal with transient 

signal properties that change over time. 

"Sample-source radiation" refers to magnetic flux emissions resulting from 
molecular motion of a sample, such as the rotation of a molecular dipole in a 
magnetic field. 

"Gaussian noise" means random noise having a Gaussian power 
distribution. 

[0052] "Stationary white Gaussian noise" means random Gaussian noise that has 
no predictable future components. 

7 10/9/03 
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[0053) 
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[0059] 



■Frequency-domain spectrum" refers to a Fourier frequency plot of a time- 
domain signal. 

"Spectral components" refer to singular or repeating qualities within a time- 
domain signal that can be measured in the frequency, amplitude, and/or phase 
domains. Spectral components will typically refer to signals present in the 
frequency domain. 

-Similar sample," with reference to a first sample, refers to the same sample 
or a sample having substantially the same sample components as the first sample. 

"Faraday cage" refers to an electromagnetic shielding configuration that 
provides an electrical path to ground for unwanted electromagnetic radiation, 
thereby quieting an electromagnetic environment. 

II. Apparatus 

Described in detail below is a system and method for detecting, processing, 
and presenting low frequency electromagnetic emissions or signals ot a sample of 

. ,, : * - L~™.m ufhi+o nr ftaiiRsian noise siqnal is 

interest, in one emoouiiMciu, a m. w » 

introduced to the sample. The Gaussian noise is configured to permit the 
electromagnetic emissions from the sample to be sufficiently detected by a signal 
detection system. Sets of detected signals are processed together to ensure 
repeatability and statistical relevance. The resulting emission pattern or spectrum 
can be displayed, stored, and/or identified as a particular substance. 

The following description provides specific details for a thorough 
understanding of, and enabling description for, embodiments of the invention. 
However, one skilled in the art will understand that the invention may be practiced 
without these details. In other instances, well-known structures and functions 
have not been shown or described in detail to avoid unnecessarily obscuring the 
description of embodiments of the invention. 

As explained in detail below, embodiments of the present invention are 
directed to providing an apparatus and method for the repeatable detection and 
recording of low-threshold molecular electromagnetic signals. A magnet.cally 
shielded faraday cage shields the sample material and detection apparatus from 
extraneous electromagnetic signals. Within the magnetically shielded faraday 
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cage, a cofl iniecte white or Gaussian nofee. a nonfenous bay holds the samp a, 
and a gradiometer detecto low*reshold moleoular eleotromagnetio signals. The 
apparatus farther factories a superconducttng quantum interferenc* de*e 
("SQUID") and a preamplifier. 

The apparatus is used by plaoing a sample within the magnetieally sh,elded 
faraday oage in close proximo to the notes coil and gradiometer. White noise © 
injected through me notee CO. and modulated unfit .he molecular electtomagne* 
signal is enhanced through stocha* resonance. The enhanced molecular 
alecttomagneuc sfanal, shietded from external interference by the faraday cage 
and the field generated by the notee at. is .ben detected and measured by the 
gra diome te r and SQUID. The signal is .hen amplified and transmuted to any 
appropriate recording or measuring equipment. 

Referring to Figure 1. mere is shown a shielding sttoctore 10 wh,ch 

.... J: o ^AnHnrtiufi wire caae 16 which is a 

includes, in an outer 10 umoi u..ov.» u .., 

m *,ne»ic sbieid and inner conducive wire cages 18 and 20 wh,ch prov,de 
e ,ec*omagne«c shielding. In another embodiment, the outer magnebc shieto . 
formed of a solid alumfaum plate material having an aluminunvnicKel altoy coafing 
and the electromagnet shielding la previded by two inner wall sbodures, each 
formed of solid aluminum. 

Referring to Figure 2, the fareday cage 10 is open a. tee top. and includes 
side openings 12 and 14. The faraday cage 10 b further comprised of mree 
copper mesh cages 16. 18 and 20. nestted in one another. Each of the copper 
mash cages 16. 18 and 20 is elecfacall, teolated from the other cages by drelectnc 

barriers (not shown) between each cage. 

Side openings 12 and 14 further comprise attenuation tubes 22 and 24 to 
provide accasa to «te interior of ttte faraday cage 10 while teolaling me intone, of 
ft. cage mom externa, sources of interference. Referring to Ffcure 3. attenuate 
tube 24 is comprised of finree copper mesh tubes 26, 26 and 30, nestled in one 
another. The exterior copper mash cages 16, 18 and 20 are each elednca ly 
exacted to one of fine copper mesh fabes 26, 28 and 30, respecWy. 
Attenuate fabe 24 is further capped with cap 32, with me cap havng hole 34. 
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Attenuation tube 22 is similariy combed of oopper mesh tubes 26. 28 and 30. 
but does not include cap 32. 

Referring again to Figure 2, a lowdensir, nonferrous sample fray 50 
mou „,ed in me interior of the farada, cage 10. The sample bay 50 H ~M so 
M « may be removed from the faraday cage 10 torough the attenuator, t*e 22 
and side opening 12. Th*e rods 52. each of wfrtoh * greater in lengtt, bran me 

o, toe attenuation tobe 22, are attach to bte sample bay 50. The brree * 
are adapted to conform to .he interior curve of the attenuation lube 22, so that the 
sample .ray 50 may be po*ned in «he center of the faraday cage 10 by resbng 
toe 1 in bre attenuabon tobe. In me Abated embodiment me sample fray 50 
and rcds 52 are made o, glass fiber epoxy. I. - * readily apparent to .hose 
MM in me art ma, me sample fray 50 and rods 52 may be made of otoe 

.ntoH in the faradav cage iu Dy 
nonferrous materials, ana we ««y 

MVxar means, such as by a single rod. 

"Referring again to Figure 2, mounted wWn me faraday cage 10 and above 
toe sample fray 50 is a cryogenic dewar 100. In the disclosed embodrmen. b» 
aewar 100 is adapted to ft within me opening a. toe top of faraday cage 10 an rs 
a Model BMD-6 Liquid Helium Dewar manufactured by Triston Technology Inc. 
The dewar 100 is constructed of a rubber epoxy composite. A gradiometer 
, ,0 w«h a very narrow beld of view b rented within me dewar 100 in pos*on so 
toa, «s freld of view encompasses me sample tray 50. In me illusfrafr* 
embodiment, toe gradiometer 110 is a brs, order axia, detechon coil. nommaM 
cenbmeter in dfcmeter. wbh a 2 % balance, and is formed from a supercondu^ 
The geometer can be any form of gradiometer excluding a planar gradrom^ 
The gradtometer 110 is connected to me input ooi, of one low ***** 
J* supemonducbng guantum interterence device W 
dfectoS ed embodiment, the SQUID is a Mode. LSQ/20 US dc SQWO 
manufactored b, Triston Technology, Inc. It will be recognized by .hose «M 
i„ toe art ma. high temperatore or aHernating curren, SQUtDs can be used wdhou. 
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«kW * ™ an T "*T ,221 dewar 100 and transfers 

„ ^ fron. 1. SQUID 1 0 o J ux ^ ^ 

« to ft. feraday cage 10 and dew r by 
the disclosed embodiment is an iFL-301-L Flux L 
Tristan Loop , 40 further amplifies and 

„ Referring to Figure . * * ^n^o***^** 
outpulstesi^^Nedfr^OreSQUtDiao^ ^ ^ ^ ^ 

— « * 8 £J r^pltcconnecsn, cabfe 144 * SQUID 
iMtotheSQUIDcontrolterloO.Thettoe p ^ controller 150 

is ^ extern* fe Ore magnehc ^^, Dcortrater150t o«» 
^^em^rscor^^als^lKeSQ ^ 

Flux Locked Loop 140, further reduorg 

Terence * *e signa. fe be measured .1 - J ^ 
beu8 edw«depan1ngfe».*»sp 1 n.and^rt I 

converters 152. a standard M ■» «• me BNC 

— rj= m — — 
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Helnd* «-*"" 60 « * ^ a „ d „an. 

„ o. 8.8 megahertt. ft «- , ^ « by 4 

«, 64 are ,«* — - «" * ^ d but shouW be shape* ami 
inch es*de. jTJI „ ^ed — me 

Hd priced by the Helmhote col B** ffi ee and 68 

^ 68 are «* arched to ™ „ „ frames pen* 

, re ,*n * *e io*er p*n - « * ^ „ » vary 

.*Mrt of the coil 62 and ^ ?0 ^ w or 

tow*****"***???! ^ ^ ^ ^ ^ 
„ epoxied omo to WM - ■- -~ ^ ^ d ^ sbe[ epoxy. 
emW «ment. the frames 66 and 68 and 9 ^ ^ ^ 

ertodinent in 

^ Refem^ to W» **■"'* ^ „ , verti cal onertafcn and an 
«, Helmh* coil *dings 62 and 6 ^ ^ ^ „ „ 

agonal no* 0 300 is ^ JJ^JL th e venial ««» « 

an0 64 c* H^^^'^^entothebo^o.^ycaga 10. 
SOOa^uss.betanV*."^ ^ onolseco „ 00from an 

MM «ed pa,r * (not «-0 ^ „ supply «, 

r^^^^-^^^r^^viaan 
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«*. o- * a balanced splitter ton an output 
MM. uoise output connect, or v ^ ^ ^ ^ ^ 

connection to the rrtsa generator. Alton ^ ^ „ 

^ C 300 would be *«* - ™ series 0( « value RF 

wnichmanyareavaitabte— y.orvraasu, 

attenuationfilters. frames suppor6n g fte coils of 

, Refers to Figure 6, a deta.1 ° ^ ^ ( fe „ 

Hrtft* transient 60 ma, be seen, toe ^ 
^.^--•Jj,, h . sub-a., 
« « disposed to *ow toe coi, - *^ „ and 6ff illustrat e toe 

^ po*n and par* to one ^ ^^jo^aa^. 

relationship v* one anotoer. ^ ^ ^ ^ ge „ erator 

m «*■"» aa * to T,TZr 40 and fe electrteally connected to toe 
SO . e*«nal to magna* *,d^e4 ^ ^ ^ „ 

He,-U transtorme, 60 - — „ lube 24 . ^ 
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[00901 The detection unit 702 is configured to handle solid, liquid, or gas samples. 
Various sample holders may be utilized in the detection unit 702. For example, 
depending on the size of the sample, a larger sample holder may be utilized. As 
another example, when the sample is reactive to air, the sample holder can be 
configured to encapsulate or form an airtight seal around the sample. In still 
another example, when the sample is in a gaseous state, the sample can be 
introduced inside the sample chamber 706 without the sample holder 720. For 
such samples, the sample chamber 706 is held at a vacuum. A vacuum seal 721 
at the top of the sample chamber 706 aids in maintaining a vacuum and/or 
accommodating the sample holder 720. 
[0091] A sense coil 722 and a sense coil 724, also referred to as detection coils, 
are provided above and below the sample holder 720, respectively. The coil 
windings of the sense coils 722, 724 are configured to operate in the direct current 

i : ^i.. cn wnhoi+y ( M7\ ranne. with a center frequency of 25 kHz 

and a self-resonant frequency of 8.8 MHz. The sense coils 722, 724 are in the 
second derivative form and are configured to achieve approximately 100% 
coupling. In one embodiment, the coils 722, 724 are generally rectangular in 
shape and are held in place by 610 fasteners. The coils 722, 724 function as a 
second derivative gradiometer. 
[0092] Helmholtz coils 726 and 728 may be vertically positioned between the 
shield 712 and the metal cage 708, as explained herein. Each of the coils 726 and 
728 may be raised or lowered independently of each other. The coils 726 and 
728, also referred to as a white or Gaussian noise generation coils, are at room or 
ambient temperature. The noise generated by the coils 726, 728 is approximately 
0.10 Gauss. 

[0093J The degree of coupling between the emissions from the sample and the 
coils 722, 724 may be changed by repositioning the sample holder 720 relative to 
the coils 722, 724, or by repositioning one or both of the coils 726, 728 relative to 
the sample holder 720. 
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(e g V or noise signal) to the coils 726, 728 and a second signal (e.g.. V or 
noise cancelled signal) to a noise cancellation coil of a superconducUng 
quantum interference device (SQUID) 206. The ampler 202 is configured to lock 
wfthout an externa, reference and may be a Peddns Elmer model 7265 DSP lock- 
in amplifier. This amplifier wort* in a "virtual mode.- where it locks to an mitral 
reference frequency, and then removes the reference frequency to allow It to run 

freely and lock to "noise." 
,0098] An analog noise generator 200 is electrically coupled to the amplifier 202. 

The generator 200 is configured to generate or induce an analog white Gauss.an 

noise at the coils 726, 728 via the amplifier 202. As an example, the generator 

200 may be a model 1380 manufactured by General Radio. 
[0099] An impedance transformer 204 Is electrically coupled between the SQU.D 

206 and the amplifier 202. The impedance transformer 204 is configured to 

. . ... _ cni nn 9nfi and amDlifier 202. 

provide impeaance maiu ■ iy 

The noise cancellation feature of the SQUID 206 can be turned on or off. 
When the noise cancellation feature is turned on, the SQUID 206 is capab.e or 
canceling or nullifying the injected noise component from the detected emiss,ons. 
To provide the noise cancellation, the first signal to the coils 726, 728 is a noise 
signal at 20 dB or 35 dB above the molecular electromagnetic emissions sought to 
be detected. At this level, the injected noise takes on the characteristics of the 
molecular electromagnetic signal through stochastic resonance. The second 
signal to the SQUID 206 is a noise cancellation signal and is inverted from the first 
signal at an amplitude sufficient to null the noise at the SQUID output (e.g., 180 
degrees out of phase with respect to the first signal). 
f00101] The SQUID 206 is a low temperature direct element SQUID. As an 
example, the SQUID 206 may be a mode. LSQ/20 LTS dC SQU.D manufactured 
by Tristan Technologies, Inc. Alternatively, a high temperature or alternating 
current SQUID can be used. The coils 722, 724 (e.g.. gradiometer) and the 
SQUID 206 (collectively referred to as the SQUID/gradiometer detector assembly) 
combined has a magnetic field measuring sensitive of approximately 5 
miccoTesla/^z. The induced voltage in the coils 722, 724 is detected and 
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Rpfpmnr.f> Subtraction: 

pong Reference Signal Subtraction (baseline noise) is performed on Y ax.s 

(amplitude) at each point (channel) along the X (time) axis. Negative Y values are 
then zeroed. 

nrnss-Correlation: 

[00116] The function calculates the cross correlation function using summation and 
integration. Since the signal is transient, the correlation function is calculated 
using direct multiplication and integration. All of the values required for the 
calculation which lie outside the source channels (data series) are taken to be 0. 
The points for which t<0 are also calculated. 

Fnnripr Significan ce Levels: 

P01171 Monte Carlo data is fitted to parametric models. Where data size N is the 

... -r. L i-/-* on noramotrir. mortals are used. For a 

only factor, univariate lauievuiv* ^ r - 

seamented FFT where segment size and overlap are additional influences, 
trivariateChebyshev polynomials are implemented. These are options selected 
under Autosignal. One could have data sets that analyze individually, or could be 
analyzed in an overlapping fashion where data set one would be analyzed, then 
the second half of data set one and the first half of data set two, then data set two, 
then the second half . 

potiq A flow diagram of the signal detection and processing performed by the 
system 100 is shown in Figure 10. When a sample is of interest, at least four 
signal detections or data runs are performed: a first data run at a time t, without 
the sample, a second data run at a time t 2 with the sample, a third data run at a 
time t 3 Wrth the sample, and a fourth data run at a time t, without the sample. 
Performing and collecting data sets from more than one data run increases 
accuracy of the final (e.g., correlated) data set. In the four data runs, the 
parameters and conditions of the system 100 are held constant (e.g., temperature, 
amount of amplification, position of the coils, the noise signal, etc.). 

[001191 At a block 300. the appropriate sample (or if it's a first or fourth data run, no 
sample), is placed in the system 100. A given sample, without injected no.se, 
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emits electromagnetic emissions in the DC-50 kHz range at an amplitude equal to 
or less than approximately 0.001 microTesla. To capture such low emiss.ons, a 
white Gaussian noise is injected at a block 301. 
P0120] At a block 302, the coils 722, 724 detect the induced voltage representative 
of the sample's emission and the injected noise. The induced voltage comprises a 
continuous stream of voltage values (amplitude and phase) as a function of time 
for the duration of a data run. A data run can be 2-20 minutes in length and 
hence, the data set corresponding to the data run comprises 2-20 minutes of 
voltage values as a function of time. 
[001211 At a block 304, the injected noise is cancelled as the induced voltage is 
being detected. This block is omitted when the noise cancellation feature of the 
SQUID 206 is turned off. 
pi* At a block 306, the voltage values of the data set are amplified by 20-50 dB, 

n-««« nmirreH at the block 304. And at a 

depending on wneinei "uo" — — 

block 308. the amplified data set undergoes analog to digital (A/D) conversion and 
is stored in the recorder 212. A digitized data set can comprise millions of rows ot 
data. 

(00123] After the acquired data set is stored, at a block 310 a check is performed to 
see whether at least four data runs for the sample have occurred (e.g., have 
acquired at least four data sets). If four data sets for a given sample have been 
obtained, then lowpass filtering occurs at a block 312. Otherwise, the next data 
run is initiated (return to the block 300). 

mm After lowpass filtering (block 312) and bandpass filtering (at a block 314) 
the digitized data sets, the data sets are converted to the frequency domain at a 
Fourier transform block 316. 

P125] Next, at a block 318, like data sets are correlated with each other at each 
data point. For example, the first data set corresponding to the first data run (e.g., 
a baseline or ambient noise data run) and the fourth data set corresponding to the 
fourth data run (e.g.. another noise data run) are correlated to each other. IF the 
amplitude value of the first data set at a given frequency is the same as the 
amplitude value of the fourth data set at that given frequency, then the correlation 
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value or number for that given frequency would be 1 .0. Alternatively, the range of 
correlation values may be set at between 0-100. Such correlation or comparison 
also occurs for the second and third data runs (e.g., the sample data runs). 
Because the acquired data sets are stored, they can be accessed at a later time 
as the remaining data runs are completed. 
I00 126] When the SQUID 206 provides no noise cancellation, then predetermined 
threshold levels are applied to each correlated data set to eliminate statistically 
irrelevant correlation values. A variety of threshold values may be used, 
depending on the length of the data runs (the longer the data runs, greater the 
accuracy of the acquired data) and the likely similarity of the sample's actual 
emission spectrum to other types of samples. In addition to the threshold levels, 
the correlations are averaged. Use of thresholds and averaging correlation results 
in the injected noise component becoming very small in the resulting correlated 
data set. 

,oo 127! If noise cancellation is provided at the SQUID 206, then the use of 
thresholds and averaging correlations are not necessary. 

I00 128] Once the two sample data sets have been refined to a correlated sample 
data set and the two noise data sets have been refined to a correlated noise data 
set, the correlated noise data set is subtracted from the correlated sample data 
set. The resulting data set is the final data set (e.g., a data set representative of 
the emission spectrum of the sample) (block 320). 

l0 oi29] Since there can be 8600 data points per Hz and the final data set can have 
data points for a frequency range of DC-50 kHz, the final data set can comprise 
several hundred million rows of data. Each row of data can include the frequency, 
amplitude, phase, and a correlation value. 

P0130] In Figures 11A and 11B. there are shown examples of sample emission 
spectrums. A Fourier plot 400 shown in Figure 1 1 A corresponds to a spectrum of 
a sample of saturated sodium chloride solution. A Fourier plot 500 shown in 
Figure 1 1B corresponds to a spectrum of a sample of an enzyme. 
P0131] Referring to Figure 16, another alternative embodiment to the systems 
described above will now be described as a system 1600. In general, alternatives 
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[00136] The random white noise generator, model 1381, manufactured by General 
Radio and described above, may be replaced by a programmable Gaussian white 
noise generator manufactured by Noise/Com. Such a generator employs two 
outputs, one inverted from the other. One output may be connected to the 
Helmholtz coil, with the other (inverted) output connected to the SQUID noise 
cancellation coil noted above. 
(00137] Likewise, as shown in Figure 19, the Tektronix digital oscilliscope noted 
above may be replaced by a two-channeled dynamic signal analyzer 1902, model 
SR 785, manufactured by Stanford Research Systems. Such a signal analyzer 
may process incoming signals by sampling multiple time domain signals and 
averaging them across multiple frequency domain FFT's. This may result in a full 
spectrum frequency domain record of all non-random signal components. Other 
changes that may be made include replacing the digital audio tape storage system 
... /n\/n\ rorvM-Hw -mru Further, a data acquisition board 

Wltfl a Uiynai woioouic uwi> v" * "i • — 

1906 manufactured by Keithley, model 3801, may be used, which works with 
software for generating histograms, as described below. 

[00138] In the alternative embodiment shown in Figure 19, a noise cancellation coil 
1908 is connected between the gradiometer and SQUID. (While a first derivative 
gradiometer is shown, a second derivative gradiometer, such as that shown in 
Figure 16, may be used.) While not shown in Figure 19, an inverted noise channel 
(inverted with respect to noise applied to the Helmholtz coils) may be applied to 
the noise cancellation coil 1908 (and may first pass through an impedance 
transformer that attenuates the noise signal by, for example, 45 dB). In an 
alternative embodiment, not shown, the noise cancellation coil may be positioned 
within the SQUID 120, between the SQUID input and output coils. 
HI. Histog ram Method nf Generating Spectraj Information 

100139] Figure 20 is a high level data flow diagram in the histogram method for 
generating spectral information. Data acquired from the SQUID (box 2002) or 
stored data (box 2004) is saved as 16 bit WAV data (box 2006), and converted 
into double-precision floating point data (box 2008). The converted data may be 
saved (box 2010) or displayed as a raw waveform (box 2012). The converted data 
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is then passed to the algorithm described below with respect to Figure 21, and 
indicated by the box 2014 labeled Fourier Analysis. The histogram can be 
displayed at 201 6. 

[00 i40i With reference to Figure 21 , the general flow of the histogram algorithm .s 
to take a discrete sampled time-domain signal and use Fourier analysis to convert 
it to a frequency domain spectrum for further analysis. The time-domain s.gnals 
are acquired from an ADC (analog/digital converter) and stored in the buffer 
indicated at 2102. This sample is SampleDuration seconds long, and is sampled 
at SampleRate samples per second, thus providing SampleCount 
{SampleDuration * SampleRate) samples. The FrequencyRange that can be 
recovered from the signal is defined as half the SampleRate, as defined by 
Nyquist. Thus, if a time-series signal is sampled at 10,000 samples per second, 
the FrequencyRange will be 0 Hz to 5 kHz. One Fourier algorfthm that may be 

_ ~ .. r-... r „.^, Troncfnrm /rfftv which has a selectable 
used is a Kaaix c rwai raoi • •• * 

freouencv domain resoMon (FFTSize) of powers of two up to 2«. An FFTSte of 
8192 is selected, to provide provides enough resolution to have at least one 
spectrum bin per Hertz as long as the FrequencyRange stays at or below 8 kHz. 
The SampteDureSPn should be long enough such that SampfeCount >(2") 
FFTSize * 10 to ensure reliable results. 

,001411 Since this FFT can only act on FFTSize samples at a time, the program 
must perform the FFT on the samples sequentially and average the results 
together to get the final spectrum. If one chooses to skip FFTSize samples for 
each FFT, a statistical error of 1 / FFTSize * 0.5 is introduced. If, however, one 
chooses to overlap the FFT input by half the FFTSize, this error is reduced to 1 / 
(081 * 2 * FFTSize) A 0.5. This reduces the error from 0.0110485435 to 
00086805556. Additional information about errors and correlation analyses in 
general, consult Bendat & Piersol, "Engineering Applications of Correlate and 
Spectral Analysis", 1993. 

P0142, Prior to performing the FFT on a given window, a data tapering filter may be 
applied to avoid spectral leakage due to sampling aliasing. This filter can be 
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chosen tan among Rectangular (no filter). Hamming, Hanning. MM 
Blackman and Blackman/Harris, as examples. 
„„, in an exemplary method, and as shown in box 2104, we have chosen 8192 
fo, the variable FFTS*e, which will be the number of time-domain samples we 
operate on a. a time, as weU as me number of discrete frequencies output by the 
FFT Note that FFTSize =8192 is the resolution, or number of bins n the range 
wWch is dictated by the sampling rate. The variable n, which dictates how many 
discrate RFFTs (Real FFTs) perforated, is set by dividing the SampleCount by 
FFTSize • 2, the number of FFT bins. In order for the algorithm to generate 
sensible results, this number n shoukt be a. least 10 to 20 (alfrtough other valves 
are possWe), where more may be preferred to pick up weaker signals. This 
tmpiies M for a given SampleRate and FFTSize, the SampleDuration must be 
long enough. A counter m, which counts from 0 to n, is initialized to zero, also as 

shown in box 2104. 

The oroaram first estabfches three buffers: buffer 2108 for FFTSize 
histogram bins, that will accumulate counts at each bin frequency; buffer 2110 for 
average power at each bin frequency, and a buffer 2112 curtaining Ihe FFTSize 

copied samples for each m. 
IOOM5 , The pmgram initializes the histograms and arrays (box 21 13) and cop.es 
FFTSize samples of the wave data into buffer 2112. at 2114, and perforata an 
RFFT on the wave data (box 2115). The FFT is norraafeed so that the highest 
amplKude is 1 (box 2116) and the average power for all FFTSize bats .s 
deterarined from the normalized signal (box 2117). For each bin frequency, the 
normalized value from the FFT at that frequency is added to each bin in buffer 
2108 (box 2118). 

mm in box 21 19 the pragram then looks at the power at each bin frequency, 
relate to the average power ovulated from above. If the power is within a 
certain (actor eosilon (between 0 and 1) of the average power, then it is counted 
and the corresponding bin is incremented in the histogram buffer at 16. Otherwtse 
it is discarded. 



xr\f\iAA\ 
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mtn Note that the average power it is comparing to is tor this FFT instance only. 
An enhanced, albeit slower algorithm might «ake two passes through the date and 
compute the average over all time before setting histogram levels. The 
comparison to epsilon helps to represent a power value that is significant enough 
for a frequency bin. Or in broader terms, the equation employing epsilon helps 
answer the question, "is there a signal at this frequency at this t,me?" If the 
answer is yes. it could due be one of two things: (1) stattonary noise whrch ,s 
landing in this bin just this one time, or (2) a real low level periodic signal wh*h wtll 
occur needy every time. Thus, the histogram counts will weed out the noise h,ts, 
and enhance the low level signal hits. So, the averaging and epsilon factor allow 
one to select the smallest power level considered significant. 
mm Counter m is incremented at box 2120, and the above process b repeated 
for each n set of WAV data unSI m is equal to n (box 2121). At each cycle, the 

• -j j- j »~ tu« coexisted hin at 2118, and each 
average power tor eacn rati is auucu u. 

hin is incremented by one when the power amplKude condfcn at 21 14 

is met. 

-M When all n cycles of data have been considered, the average power ,n each 
bin is determined by dMding the total accumulated average power in each btn by 
„ the total number of cycles (box 2122) and the results displayed (box 2123). 
Except where structured noise exists. e.g„ DC = 0 or at multiples of 60 Hz, the 
average power in each bin will be some relatively low number. This is indicated m 
the plots shown at Figures 22A-D (the holograms generated at 400, 600, 700, and 
900 mV) The plots of Figures 22A-22D show only a portion of the histogram bins, 
namely a specttum tan 7953Hz through 8533Hz. As shown in Figures 22A and 
22B no stochas«c event is visible at 400 mV or 600 mV of injected noise, 
respectively. However, as shown in Figure 22C, a. 700 mV. a visible stochastto 
even, is evident. Thereafter, as shown in Ffcure 22D, a. 900 mV, the stochaste 
event is lost. 

[00150] The histogram produced by the above steps contains, in each bin, a count 
between 0 and n of the number of times that the power at that frequency was 
above (epsilon * the average power for that whole FFT output). If a bin count .s 
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incremented due to —red noise, that noise « be distnbuted across all me 
frequency bins over time, thus no. adding up to much in a given to. if there . 
consent signai a. a given frequency, it « be present a, each of the n ume slrcas 
and thus nave a bin count approaching n. Large amplitude noise, such as sorry 
nedz and its harmonics have both high to counts as weH as high average power. 
We can differentiate between these frequencies, and the ones we are interested rn 
that have low average power, but high bin counts. 

Figure 22A-22D show histograms generated by the method at four 
different noise power Inpute. As shown, the program may display average power 
at each frequency as a vertical bar. The histogram bin counts may be represented 
as a connected upper line. « the power b ordered W (e.g. less man 
average/3), and me histogram has a certain count, men a connecffng Ime may 
become observable between the peak of a power bar and a peak of a htetogram 

x:„ lii^oiu candidates for low energy 

bar. Bins highlighted Dy me wmmmhv — 

mAlon ilar sDectra. 

"""""itcan be appteciated from Figures 22A-22D and from the above, that there 
are two settings of note used m generating a meaning*, histogram, mat rs. a 
tetogram M shows stochasfc resonance effects related to a samp* bemg 
interrogated. The first is the power level of Gaussian vrhite noise supplied to (he 
aample. If this level istoo low, me noise level is no, sufficient create stocha* 
resonance andme bin histogram reflectenofceonly. If the power input is too hrg , 
toe average power level calculated tor eacb bin is high and st«hastx even* 

cannot be distinguished. 
mm The second seeing is «te value of eps»on. This value determines a power 
valuethat-beusedtodistmguish^eventoveraveragevalue. Atavalueoft, 
no everts will be detected, ste power will never be greater than average power. 
As epsikm apaches zero, virtual every value will be placed in a bin. Between 
0 and 1, and typfcally at a value mat gives a number of bin counts between about 
20-50% of total bin counts for MM r»ise, epsilon will have a rnaomum 
^gW character; meaning me stoohasho resonance events will be most h,Wy 
favored over pure noise. 

10/9/03 
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[00154) Therefore, in practicing the invention, one can systematically increase the 
power gain on the noise input, e.g.. in 100 mV increments between 0 and 1 V, and 
at each power setting, adjust epsilon until a histogram having well defined peaks is 
observed. Where, for example, the sample being processed represents a 20 
second time interval, total processing time for each different power and epsilon will 
be about 25 seconds. When a well-defined signal is observed, either the power 
setting or epsilon or both can be refined until an optimal histogram, meaning one 
with the largest number of identifiable peaks, is produced. 
100155] Under this algorithm, numerous bins may be filled and associated 
histogram rendered for low frequencies due to the general occurrence of noise 
(such as environmental noise) at the low frequencies. Thus, the system may 
simply ignore bins below a given frequency (e.g., below 1 kHz), but still render 
sufficient bin values at higher frequencies to determine unique signal signatures 
between sampies. 

rooiss, Alternatively, since a purpose of the epsilon variable is to accommodate 
different average power levels determined in each cycle, the program could itself 
automatically adjust epsilon using a predefined function relating average power 
level to an optimal value of epsilon. 

100157] Similarly, the program could compare peak heights at each power setting, 
and automatically adjust the noise power setting until optimal peak heights or 
' character is observed in the histograms. 

(ooisb] Although the value of epsilon may be a fixed value for all frequencies, it is 
also contemplated to employ a frequency-dependent value for epsilon. to adjust 
for the higher value average energies that may be observed at low frequencies, 
e g., DC to 1 ,000. A frequency-dependent epsilon factor could be determined, for 
example, by averaging a large number of low-frequency FFT regions, and 
determining a value of epsilon that "adjusts" average values to values comparable 
to those observed at higher frequencies. 

P0159, Referring to Figures 23A-23C. an example of a user interface for generating 
histograms is shown. A slider bar 2302 determines the length of a sample 
waveform segment, such as up to 300-600 seconds, and allows a user to 
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effectively scroll within a waveform. A box 2304 allows the user to set a Nyquist 
frequency, such as 5, 10 or 20 kHz, and also provided is an adjacent reset button. 
A slider bar 2306 allows the user to move the baseline for histograms, while a 
60Hz checkbox 2308 allows the user to identify the 60Hz bin and all related 60Hz 
harmonics with vertical lines (as shown in Figure 23C). When an acquire button 
2312 is selected, the software generates or acquires a waveform from a sample, 
such as that shown in Figure 23B. When an fit button 2310 is selected, the 
software generates a histogram plot, such as that shown in Figure 23C. 

IV. Methods and Applications 

poieoi This section describes the use of the apparatus described above for 
interrogating a sample, and a variety of applications of the apparatus in 
characterizing a sample and in detecting sample components. Also disclosed, in 

cinnaturfi nr data 

accordance with the invention, is a iuw-ne M u C ..^ o K - r- — 

set by which a sample can be characterized, and a time-domain signal of a 
sample, used, for example, in generating the sample spectroscopic signature. 

A. MRthnd of inte rrogating a sample 

[00161] An objective of the method of the invention is to generate spectroscopic 
information relating a sample being interrogated. As will be seen, the information 
may be in the form of a spectral plot, in a selected low-frequency spectral range, 
or a data set which identifies low-frequency spectral components characterizing 
the sample, or actual identification of a sample or sample components, based on 
the characteristic frequencies identified for the sample. 

[00162] The sample may be any material having atomic or molecular components, 
e.g., ionic salt components or molecular compound in ionized or nonionized form, 
or protonated or non-protonated form, that has molecular rotation, and preferably a 
dipole moment such that molecular rotation in a magnetic field, e.g., the earth's 
magnetic field, is effective to produce a low-frequency electromagnetic emission. 
The sample is typically a liquid sample, but may be gaseous or solid or semi-solid 
as well, as long as at least one component of the sample has one or more 
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rotation* degrees of fteedont. Typical samples are aqueous ot organic solutions 
having one ot mote solute component, which may be the sample matenal of 

interest, dissolved in the solvent. 

The sample is placed in a suitable vessel, ptefetabty one such as Pytex 
glass that has little obsetvable low-frequency spectral components, and the vessel 
is then positioned in the apparatos container as descdbed in Secta, II. Wh the 
sample posifoned in the apparatus container, the Gaussian noise generator ,s 
activated to inject Gaussian noise into the sample. The amptade (mean 
amplitude) of Gaussian noise injected is preferably sufficient to produce non- 
stationary composite ftne-domain signal components. This may be done, for 
example, using an oscilloscope with a Fourier transform capabilfty, and obsennng 
the frequency-domain signal in a suitable range, e.g., 20MOOHZ window. A 
suitable noise level is selected when detectable frequency components are first 
observed. 

nurino noise injection, the recording device records a bme-domatn 

electromagnetic sfenal from the detector over a preset time intetval. The recording 
intetval may be relatively short, e.g., 3M0 seconds, or may be several minutes or 
m ore, depending on the M spectral resoMon required. The signals recorded 
ara stored in a suteble sfcnal storage device, e.g., a tape or hard disc. fo. use m 
later signal processing operations now to be described. 

mm in general, it is desirable to enhance sample signal components by cross- 
coveting the sample timeOomain signal recorded with a second tota> 
signal of the seme sample or, less preferably, an identical sample or a sample 
having the same sample components of interest The receding time for the 
second sfcnal is preferably the same as for the M sfcnal. The two signals are 
cosseted using a standard cross^rrelaton algorithm in the time domaa 
This results in a spreadsheet or spectrom idenWying the signal spectral 
components that are common in bo* signals that hold up over time, and a 
correlate value for each component wh» measures the relationship botween 
spectral components common to both signals. 
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l0 oi66i The improvement in spectral resolution obtained by the signal cross- 
correlation is seen in the Figures 12A and 12B, and Figures 13A and 13B. The 
Figures are Fast Fourier transforms of a first time-dependent signal in the 
frequency domain (Figures 12A and 13A) or a fast Fourier transform of the first 
and second cross-correlated frequency-domain spectra (the spreadsheet referred 
to above) (Figures 13B and 13B) to plot the spectral components in the frequency 
domain, and in the spectral range of 500-530 Hz. 

[00167, Comparing Figures 12A and 12B for a NaCl sample, it is seen that the 
cross-correlation signal processing significant enhances signal-to-no.se ratio, 
brings out much more detail in the sample-specific spectral component whose 
peak is at 522.5 Hz, and also produces a signfficantly refined peak position. 
Similar (and exemplary) results were observed for the alkyl ether sulfate sample, 
whose spectral features in the 500-530 Hz range are seen in Figures 13A and 13B 

j - : io ^crw-tiv/oiv As with the NaCl sample, the 

for unoorreiaiea auu wnewou oi a ..«— , r ~ *■ 

soectrum derived from the correlated signals gave much lower signal-to-noise 
ratio much more detail and information as to sample-specfic spectral 
components. The signal correlation can also be applied, conventionally, to 
produce a spreadsheet relating frequency and phase (rather than frequency and 
amplitude). 

[00168] The correlated time-domain spectrum from above is plotted in the frequency 
domain by applying a fast Fourier transform to the spectrum, where the spectral 
correlations values are represented in the y axis as amplitudes. The plot is wrth.n 
the frequency range DC to 50 kHz, preferably in the region DC to 6.5 kHz. As will 
be seen below, the dominant spectral features of many samples are found m the 
100-1.500 Hz range, particularly the 500-550 Hz range; and thus the spectrum 
generated may be confined accordingly, e.g.. in the 500-530 Hz range. The FTT 
is carried out by a well known FTT algorithm. The correlated time-domain s.gnal 
may also, or alternatively, be transformed to phase-domain or amplitude or 
magnitude domain signals, to extract signal information related to phase or 
amplitude components of the sample spectrum. 
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(00169] Once the frequency-domain spectrum is generated, either by the cross- 
correlating or FFT step, the spectrum is used to identify one or more low- 
frequency signal components that are characteristic of the sample be.ng 
interrogated. This step may be performed by the user from direct viewing or by 
computer analysis of the spectra. 

Figures 14A-14F show spectral features for the samples deionized water 
(Figure 14A), saturated NaCI (Figure 14B). 1% NaCI in deionized water (Figure 
14C) saturated NaBr (Figure 14D), alkyl ether sulfate (Figure 14E), and empty 
sample vessel (Figure 14F), all in the spectral range of about 500-530. As seen, 
each sample has distinctive spectral components characterized by one or more 
peaks at well defined frequencies. 

B r.haracterizina a sample 

Accordance with another aspect of the invention, the method above is used 
' to generate a data set of tow-frequency spectral components of a given samp®, 

nf tho material. 

also referred to herein as a iow-Trequeiiuy & .y o - a 

[001721 The 500-530 spectral range shown for the several samples above was 
selected to illustrate having prominent spectral features in the various samples. In 
order to obtain a more complete data set of spectral frequency components of a 
sample, spectral components over a wider frequency range, e.g., 100-1,500 Hz, 
should be determined. In one aspect, the invention includes a data set of spectral 
components associated with a given sample material, e.g.. a solvent, gas, or 
solute component of a solution. The data set includes a list of the low-frequency 
spectral components of the sample, e.g.. in the 100-1,500 Hz range, whose cross- 
spectral correlates have a se.ected statistical measure above background 
spectral noise, or selected ones of these components that are unique to the 
sample. 

ipo,,,, A variety of signal-analysis methods may be employed in generator* the 
low-frequency data set for a given sample. In one exemplary method, a cross- 
correlated sample signal spectrum is compared with a cross^elated notse (no 
sample) signal. The algorithm next advances incremental, e.g.. in 0.1 Hz intetvals 
aorcss the cnossWated sanple spectrum and the cross-corrected no«e 
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spectrum, looking at the correlation value at each frequency point, and subtracts 
the noise correlation from the sample correlation at that point, to yield a frequency 
plot of corrected correlation values. These values will be relative to a particular 
sample, and depend, for example, on the relative amplitude of any noise 
component. 

mm In general, frequency components having a higher correlation value 
(relative to the other values in that sample) will tend to hold up (be observed) over 
many interrogations of the same sample. To identify those that do hold up, the 
frequency components observed for the sample over two or more sample sets, 
each obtained as above, are compared, and only those that are seen in two (or 
more, if available) sets are taken as valid components of the data set for that 
sample. In the tables below, data sets for several samples (as identified in the 
tables) are given along with the correlations determined from a single sample 

-i. ■^A-.^taA in itaiirs rtvnicallv having, the smaller 

interrogation, hwm vaiuco — — «■>. 

correlation values) were found not to hold up in multiple data sets from the same 
sample material. 

pom] Thus, for example, for the saturated NaCI sample in Table 1, spectral 
components at 522.58, 523.12, 523.47, and 523.85 Hz correlate from sample to 
sample, and would form a data set for the sample in the frequency range 500-530 
Hz. Additional members of the data set may be included In an expanded 
frequency range. 

,00176, Similarly, for the amino acid sample of Table 3, the data set would include 
components at 262.93, 257.81, 257.23, 536.68, 448.05, 531.37, 528.80, 593.44, 
588.68, 583.74, 578.61, 769.59. and 744.14 in the frequency range of between 
about 250 and 1,400 Hz. The greater spectral composition of the amino acid 
sample, relative to NaCI, presumably reflects in part, the greater complexity of the 
sample molecule. 
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Spectra 
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24 


744.14 


.039 ^ 



The data above demonstrates that both simple and more complex 
molecular samples can be characterized in terms of unique low-frequency spectral 
components. The data set associated with a given sample material may also 
include (as shown in the tables) the associated correlations values of the spectral 
components. The data set may be used for example, in identifying components in 
an unknown sample and/or for estimating the relative concentrations of a material 
in a sample. The use of the method for identifying low-concentration components 
in a sample is discussed in the next section. 
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C. Identif ying co ™P nnRnts in a sample 
[00178! It is often desirable to detect sample components, e.g., trace contaminants, 
present in a multi-component sample material, such as a liquid sample with 
unknown contaminants, or other samples capable of holding or supporting a 
contaminant that it is desired to detect. 
[00179] An analytical method for detecting a component of a sample, in accordance 
with another aspect of the invention, includes first identifying the low-frequency 
sample spectral components of a sample (i) in a selected frequency range 
between DC and 50 kHz (ii) whose cross-spectral correlations have a selected 
statistical measure above background spectral noise, as described above. 
[00180] The sample spectral components are then compared with characteristic 
low-frequency spectral components of known compounds suspected of being 
present in the sample. In a typical example, the sample components are 

. . : _ -x «... n f o^h Af thp rnmnonents suspected of being in 
compareu ayamoi mo «ou. — 

the sample and which one desires to detect. A components, e.g., compound is 
identified as being present in the sample if one or more of its characteristic low- 
frequency spectral components correspond to one or more low-frequency spectral 
components of a known sample. 

I00181] As shown in the set of Figures 11A-11F, detection of a compound (an 
amino acid) can occur at very low levels, e.g.. in the parts per billion range or 
lower. In particular, even at a dilution of 1:10 billion w/v, a characteristic spectral 
component at about 531 Hz is observed. The Figures demonstrate that signal 
amplitude, corresponding spectral component correlation, does decline with 
increasing compound dilution. However, the loss in signal amplitude at low 
concentration can be compensated for by extending the recording time, in this 
example, from 50 seconds for the first group of Figures to 4.25 minutes for the 
most dilute sample (Figure 1 1F). 

[00182] Where, as in the above example, the spectral component amplitude 
declines with decreasing concentration, the amount of compound can be 
estimated on the basis of signal amplitude, assuming that the data set for the 
compound also includes concentration dependent amplitude information. 

, 0 _ 10/9/03 
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[00183] It has also been observed in some cases that the frequency of the 
characteristic spectral components may shift by as much as 3 Hz in a systematic 
fashion with changes in concentration. For such compounds, the amount of 
material present in a sample can be estimated by changes in amplitude and/or 
frequency shift in one or more of the spectral components. It will be appreciated 
for materials showing a concentration dependent frequency shift that a data set for 
that compound could include concentration-dependent frequencies as well as 
concentration-dependent amplitudes for particular components. 

D. Timfi-domain signals 
[00184] In still another aspect, the invention includes a time-domain signal 
associated with a material of interest. The time-domain signal, and its method of 
production, have been discussed above. Briefly, the signal is produced by placing 
the sample of interest in a container having both magnetic and electromagnetic 
shielding, injecting Gaussian noise into the sample; and recording an 

. . .. ...... j : _ „i ^mnneoH n{ camnlfi source radiation 

eiectromagneuu uiir»-uuiihuii oiynai whk«— — — «-•- -- 

superimposed on the injected Gaussian noise. 
[00185] The signal may be used to characterize a sample, much as a spectral 
component data set is used to characterize a material. Alternatively, the signal 
may be used for generating a low frequency signal signature of spectral 
components associated with a material of interest. The signal signature can be 
generated, also as described above, by (i) cross correlating the time-domain signal 
recorded with a second time domain signal separately recorded from the same or 
similar sample, to produce a frequency domain spectrum in a frequency range 
within DC to 50 kHz. 
Conclusion 

[00186] Unless the context clearly requires otherwise, throughout the description 
and the claims, the words "comprise," "comprising," and the like are to be 
construed in an inclusive sense as opposed to an exclusive or exhaustive sense; 
that is to say, in the sense of "including, but not limited to." Words in the above 
detailed description using the singular or plural number may also include the plural 
or singular number respectively. Additionally, the words "herein," "above." "below" 
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and words of similar import, when used in this application, shall refer to this 
application as a whole and not to any particular portions of this application. When 
the claims use the word "or" in reference to a list of two or more items, that word 
covers all of the following interpretations of the word: any of the items in the list, all 
of the items in the list and any combination of the items in the list. 
(00187] The above detailed descriptions of embodiments of the invention are not 
intended to be exhaustive or to limit the invention to the precise form disclosed 
above. While specific embodiments of, and examples for, the invention are 
described above for illustrative purposes, various equivalent modifications are 
possible within the scope of the invention, as those skilled in the relevant art will 
recognize. For example, while processes or steps are presented in a given order, 
alternative embodiments may perform routines having steps in a different order, 
and some steps may be deleted, moved, added, subdivided, combined, and/or 

... . - . r 11 x u« irv^iomnntoH in a varietv of different ways. 

moamea. cacn ui mwc oi^o ...o ; ~w 

Also, while these steps are shown as being performed in series, these steps may 
instead be performed in parallel, or may be performed at different times. 
,001881 The teachings of the invention provided herein can be applied to other 
systems, not necessarily the system described herein. These and other changes 
can be made to the invention in light of the detailed description. The elements and 
acts of the various embodiments described above can be combined to provide 
further embodiments. 

1001891 All of the above patents and applications and other references, including 
any that may be listed in accompanying filing papers, are incorporated herein by 
reference. Aspects of the invention can be modified, if necessary, to employ the 
systems, functions and concepts of the various references described above to 
provide yet further embodiments of the invention. 

[001901 These and other changes can be made to the invention in light of the above 
detailed description. While the above description details certain embodiments of 
the invention and describes the best mode contemplated, no matter how detailed 
the above appears in text, the invention can be practiced in many ways. Details of 
the signal acquiring and analysis system may vary considerably in its 
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implementation details, while still be encompassed by the invention disclosed 
herein. As noted above, particular terminology used when describing certain 
features or aspects of the invention should not be taken to imply that the 
terminology is being re-defined herein to be restricted to any specific 
characteristics, features or aspects of the invention with which that terminology is 
associated. In general, the terms used in the following claims should not be 
construed to limit the invention to the specific embodiments disclosed in the 
specification, unless the above Detailed Description section explicitly defines such 
terms. Accordingly, the actual scope of the invention encompasses not only the 
disclosed embodiments, but also all equivalent ways of practicing or implementing 
the invention under the claims. 

While certain aspects of the invention are presented below in certain claim 
forms, the inventors contemplate the various aspects of the invention in any 
pumhor 0 f claim forms Fnr example., while only one aspect of the invention is 
recited as embodied in a method claim format, it may likewise be embodied in a 
computer-readable medium claim format. Accordingly, the inventors reserve the 
right to add additional claims after filing the application to pursue such additional 
claim forms for other aspects of the invention. 
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CLAIMS 



I/We claim: 



1. Apparatus for interrogating a sample that exhibits low-frequency 

molecular motion, comprising 

a container adapted for receiving said sample, said container having both 

magnetic and electromagnetic shielding, 
an adjustable-power source of Gaussian noise for injection into the sample, 

with the sample in said container, 
a detector for detecting an electromagnetic time-domain signal composed 

of sample source radiation superimposed on the injected Gaussian 

noise, and 

an electronic computer aaapiea iu iweive u« u.. 

detector, and to process the signal to generate a spectral plot that 
displays, at a selected power setting of the Gaussian noise source, 
low-frequency spectral components characteristic of the sample in a 
selected frequency range between DC and 50 kHz. 

2. The apparatus of claim 1 , wherein said electronic computer includes 
a signal analyzer that functions to (i) calculate a series of Fourier spectra of the 
time-domain signal over each of a plurality of defined time periods, in a selected 
frequency range between 100 and 50 kHz, and (ii) average the Fourier spectra. 

3. The apparatus of daim 2, wherein said calculating includes 
calculating at least five Fourier spectra, each taken over a 1-5 second time-domain 
interval. 

4. The apparatus of claim 1 , wherein said electronic computer includes 

machine readable code operable to: 

(i) store a time-domain signal of the sample over a sample-duration time T; 

aa 10/9/03 
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(ii) select a sampling rate F for sampling the time domain signal, where F* T 

is the total sample count S, F is approximately twice the frequency 
domain resolution f of a Real Fast Fourier Transform of the time- 
domain signal sampled at sampling rate F, and S>f * n, where n is at 
least 10, 

(iii) select S/n samples from the stored time domain signal and performing a 

Real Fast Fourier Transform (RFFT) on the samples, 

(iv) normalize the RFFT signal and calculating an average power for the 

signal, 

(v) place an event count in each of f selected-frequency event bins where 

the measured power at the corresponding selected frequency > 
average power * e obtains, where 0< £<1, and is chosen such that 
the total number of counts placed in an event bin is between about 
20-50% nf thp maximum possible bin counts in that bin, 

(vi) repeat steps (iii-v) times, and 

(viii) generate a histogram that shows, for each event bin f over a selected 
frequency range, the number of event counts in each bin. 

5. The apparatus of claim 4, wherein said machine readable code is 
further operable to, in step (iv) place the normalized power value from the RFFT in 
f corresponding-frequency power bins, and in step (viii) (a) divide the accumulated 
values placed in each of the f power bins by n, to yield an average power in each 
bin, and (b) display on the histogram, the average power in each bin. 

6. The apparatus of claim 4, wherein said machine readable code is 
further operable to, in step (viii), identify those bins in the histogram that have an 
event count above a given threshold and a average power. 

7. The apparatus of claim 1, wherein the source of Gaussian noise 
includes an adjustable-power Gaussian noise generator and a Helmholz coil which 
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is contained within the magnetic cage and the Faraday cage, and which receives a 
selected noise output signal from the noise generator in the range 100 mV to 1 V. 

8. The apparatus of claim 7, wherein said injector is designed to inject 
Gaussian noise into the sample at a frequency between DC and 2 kHz. 

9. The apparatus of claim 1, wherein said detector is a first-derivative 
superconducting gradiometer which outputs a current signal, and a SQUID 
operatively connected to the gradiometer to convert the current signal to an 
amplified voltage signal. 

10. The apparatus of claim 9, wherein said container is an attenuation 
tube having a sample-holding region, a magnetic shielding cage surrounding said 

. - . . a ...Mi;* iha msnnai'ir chielriinn r.ane and also 

region, anu a raiauay uaye wmainou im ~ "—a 

surrounding said region, the source of Gaussian noise includes a Gaussian noise 
generator and a Helmholz coil which is contained within the magnetic cage and 
the Faraday cage, and which receives a noise output signal from the noise 
generator, and which further includes, for use in removing stationary noise 
components in the time-dependent signal, a signal inverter operatively connected 
to the said noise source and to said SQUID, for receiving Gaussian noise from the 
noise source and outputting into said SQUID, Gaussian noise in inverted form with 
respect to the Gaussian noise injected into the sample. 

11. A method for interrogating a sample that exhibits low-frequency 

molecular motion, comprising 

placing the sample in a container having both magnetic and 

electromagnetic shielding, 

(a) injecting Gaussian noise into the sample at a selected noise 
amplitude; 
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(b) recording an electromagnetic time-domain signal composed of 

sample source radiation superimposed on the injected 
Gaussian noise, 

(c) generating a spectral plot that contains, at a selected power 

setting of the Gaussian noise source, low-frequency, sample- 
dependent spectral components characteristic of the sample 
in a selected frequency range between 100 and 50 kHz, and 

(d) repeating steps (a)-(c) at different selected noise amplitudes until 

a plot showing a maximum or near maximum number of 
spectral components characteristic of the sample are 
generated. 

[ci2] 12. The method of claim 11, wherein said generating includes (i) 

onoMr* nf *ho *imp-Hnmflin aianal over each of a 

plurality of defined time periods, in a selected frequency range between 100 and 
50 kHz, and (ii) averaging the Fourier spectra. 

[ C 13] 1 3. The method of claim 1 1 , wherein said calculating includes 

(i) storing a time-domain signal of the sample over a sample-duration time 

T; 

(ii) selecting a sampling rate F for sampling the time domain signal, where 

F* T is the total sample count S, F is approximately twice the 
frequency domain resolution f of a Real Fast Fourier Transform of 
the time-domain signal sampled at sampling rate F, and S>f * n, 
where n is at least 10, 

(iii) selecting S/n samples from the stored time domain signal and 

performing a Real Fast Fourier Transform (RFFT) on the samples, 

(iv) normalizing the RFFT signal and calculating an average power for the 

signal, 

(v) placing an event count in each of f selected-frequency event bins where 

the measured power at the corresponding selected frequency > 

[/Final CIPAppln.docl 



average power * e obtains, where 0< e<1, and is chosen such that 
the total number of counts placed in an event bin is between about 
20-50% of the maximum possible bin counts in that bin, 
(vi) repeating steps (iii-v) times, and 

(viii) generating a histogram that shows, for each event bin f over a selected 
frequency range, the number of event counts in each bin. 

14. The method of claim 13, which further includes, in step (iv) placing 
the normalized power value from the RFFT in f corresponding-frequency power 
bins, and in step (viii) (a) dividing the accumulated values placed in each of the f 
power bins by n, to yield an average power in each bin, and (b) displaying on the 
histogram, the average power in each bin. 

•ic Tk„ mathnH rxf Haim 14 whinh further includes, in steD 

identifying those bins in the histogram that have an event count above a given 
threshold and a average power. 

16. A method of characterizing spectral emission features of a sample 
material, over a selected frequency range R, comprising 

(i) storing a time-domain signal of the sample over a sample-duration time 

T; 

(ii) selecting a sampling rate F for sampling the time domain signal, where 

F* T is the total sample count S, F is approximately twice the 
frequency domain resolution f of a Real Fast Fourier Transform of 
the time-domain signal sampled at sampling rate F, and S>f * n, 
where n is at least 10, 

(iii) selecting S/n samples from the stored time domain signal and 

performing a Real Fast Fourier Transform (RFFT) on the samples, 

(iv) normalizing the RFFT signal and calculating an average power for the 

signal, 
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(v) placing an event count in each of f selected-frequency event bins where 

the measured power at the corresponding selected frequency > 
average power * e obtains, where 0< e<1, and is chosen such that 
the total number of counts placed in an event bin is between about 
20-50% of the maximum possible bin counts in that bin, 

(vi) repeating steps (iii-v) times, and 

(viii) generating a histogram that shows, for each event bin f over a selected 
frequency range, the number of event counts in each bin. 

17. The method of claim 16, which further includes, in step (iv) placing 
the normalized power value from the RFFT in f corresponding-frequency power 
bins, and in step (viii) (a) dividing the accumulated values placed in each of the f 
power bins by n, to yield an average power in each bin, and (b) displaying on the 

aorh Win 

Mioiuyicun, UK? atwuyv ... ~ 

18. The method of claim 17, which further includes, in step (viii), 
identifying those bins in the histogram that have an event count above a given 
threshold and an average power. 

19. The method of claim 18, wherein R, expressed in Hz, is 
approximately equal to f, and the sample rate F, expressed in samples/second, is 
approximately 2f. 

20. The method of claim 19, wherein the method detects low-frequency 
emission events related to molecular in a sample, and wherein R includes at least 
the frequency range 1 -5 kHz. 

21. A low-frequency spectral signature associated with a material of 
interest comprising 

a list of frequency components in the DC-50 kHz frequency range that are 
generated by the method of claim 16. 

_aq 10/9/03 
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22. The spectral signature of claim 21 , wherein the frequencies in said 
list are identified from a histogram of the number of sample-dependent stochastic 
events occurring at each of a plurality of spectral increments within a selected 
frequency range between DC and 50 kHz. 
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SYSTEM AND METHOD FOR CHA RACTERIZING A SAMPLE BY LOW- 
FRFQUENCY SPECTRA 

ABSTRACT OF THE DISCLOSURE 

A method and apparatus for interrogating a sample that exhibits molecular 
rotation are disclosed. In practicing the method, the sample is placed in a 
container having both magnetic and electromagnetic shielding, and Gaussian 
noise is injected into the sample. An electromagnetic time-domain signal 
composed of sample source radiation superimposed on the injected Guassian 
noise is detected, and this signal is used to generate a spectral plot that displays, 
at a selected power setting of the Gaussian noise source, low-frequency spectral 
components characteristic of the sample in a selected frequency range between 
DC and 50 kHz. In one embodiment, the spectral plot thai is generated is a 
«f etnrhactir resnnanre fivents over the selected frequency range. 

I I Ulll VI WW i i m ■ — 

From this spectrum, one or more low-frequency signal components that are 
characteristic of the sample being interrogated are identified. 
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Example for Stochastic Generation 
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400 mV noise signal applied. No stochastic event visible. 
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600 mV noise signal applied. No stochastic event visible 
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700 mV noise signal applied. Stochastic event visible. 
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900 mV noise signal applied. No stochastic event visible. 
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